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A practical quantum-field treatment is developed for sys- 
tems endowed with an effective mass function depending on 
both space and time and a schematic application illustrates 
the quantitative importance of quantum fluctuations in the 
dynamics of disoriented chiral condensates. 



I. INTRODUCTION 

The prospect of observing manifestations of chiral sym- 
metry restoration in high-energy nuclear collisions has 
stimulated significant research activity over the past sev- 
eral years JD]. In this endeavor, the most popular theoret- 
ical tool has been the linear a model, which in its simplest 
form describes an interacting 0(4) field with isoscalar (a) 
and isovector (tt) components. Within this framework, a 
number of instructive dynamical simulations have been 
carried out at the level of classical fields for the types 
of environment expected in high-energy collisions ^-|) . 
Furthermore, quantum-field treatments have been devel- 
oped and applied, albeit within the mean-field approxi- 
mation and for uniform environments 

Relative to the classical studies, these latter treatments 
take account of the ever present quantum fluctuations 
and it appears that these may contribute significantly to 
the emerging signals ||. Therefore, it is of interest to fur- 
ther develop the quantal treatment so that more refined 
scenarios can be addressed dynamically. The present pa- 
per describes how the quantum-field treatment may be 
extended to non-uniform environments in a manner that 
allows numerical simulations for the types of scenario 
relevant to the study of disoriented chiral condensates. 
Moreover, the important role of quantum fluctuations is 
illustrated by means of a concrete application. 

The highly excited systems formed in high-energy nu- 
clear collisions expand very fast and the environments 
cool down correspondingly rapidly. As a consequence, 
the effective pion mass decreases quite quickly from its 
initially high thermal value towards its free value. Fur- 
thermore, the non-equilibrium relaxation dynamics of the 
chiral order parameter gives rise to a fairly regular oscil- 
latory modulation of the mass with a frequency near half 
the value employed for the a meson H . Since the expan- 
sion is predominantly longitudinal, it appears (3||| that 
the cooling rate is insufficient to produce a quench 0] , in 



which the square of the pion mass turns negative, causing 
the softest modes to become unstable and thus display 
an exponential population growth. While perhaps un- 
fortunate from the physics perspective, this feature does 
simplify the physical analysis and improves the prospect 
for developing test-particle transport treatments Jl(J . (A 
concrete impression of the typical time evolution of the 
effective pion mass in the bulk of the rapidly expanding 
system can be gained from Fig. [I].) 

A further complication encountered in nuclear collision 
experiments is the non-uniformity of the environment re- 
sulting from the finite extension of the systems formed. 
In the mean-field framework, this causes the mass to be 
dependent on position as well, which then couples the 
elementary modes and thus complicates the treatment 
significantly. (A typical spatial profile of the mass func- 
tion expected in nuclear collisions is shown in Fig. |[) 



II. PROPAGATION 

With such environments in mind, we now turn to the 
formal developments. In the mean-field approximation, 
the Hamiltonian density operator is given by 



W(x) = i[p(x) 2 + (Vg(x)) 2 + M 2 (x,t)q(x) 2 



(1) 



where g(x) is the field operator and p(x) is its time 
derivative. Free pions provide a natural basis but, in 
order to bring the problem onto a canonical form, each 
pair of traveling waves | ± k) is replaced by a correspond- 
ing pair of standing waves | ± k) |J . With an arbitrary 
sign convention for the quantum numbers k — ±k, the 
associated canonical eigenfunctions are 
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y/2 cos(k • x) , k, > , 

1 , k = , (2) 

— \pl sin(K • x) , k < . 



The canonical modes are orthonormal and complete, 
/dx(K|x)(x|K') = S KK ,, Ek< x I«X«I x '> = The 
canonical annihilation operator is a K , with [a K ,cc,] = 
Skk 1 , and the dispersion relation is lo 2 = k 2 + m 2 . 
The field operators can then be expanded as follows, 

g(x) = £(x|k) q K = [ fA W«a + ^ A (x)*4 

K A 

p(x) = £(x|k) p K = J2 [G A (x)a A + G A (x)*4 

A 



(4) 
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The mode functions are governed by the usual field equa- 
tion of motion, which can readily be solved numerically, 

F A (x) = G A (x), G A (x) = [A-„ 2 (x)]F A (x) . (5) 

The appropriate initial conditions are 



^=(x|A>,G A (x,t i )--^^(x|A). (6) 



Using the canonial representation of the mode functions, 
F A (x) = £<x|k) F a , G a (x) = 5>|k) G a , (7) 

we thus obtain the evolution of the canonical operators, 

(8) 



$«(*) = E [ F - (*) «a + 4 

A 



(9) 



III. ANALYSIS 

It is instructive to express the Hamiltonian in terms of 
the canonical operators p K and q K , 

H(t) = J dx W(x) = + q K M RR ,q K , . (10) 

Since the mass matrix is real and symmetric, M KK '(t) — 
(/«|/i 2 (x, t) — A\k'), it can be diagonalized at any given 
time t by a real, unitary matrix S(t), leading to 



Hit) = IE[^ + J&c& 

K 



(11) 



Here the correlated canonical field operators are given by 
Ok = J^kSkkW , Pk = X^'W*) > (I 2 ) 

K K 

and they satisfy [Qk, Pk'] = i$KK' • At the given time t, 
the corresponding correlated eigenfunctions are 



(x|K)=J>|fs) S KK (t) . 



(13) 



If the eigenvalues are positive, Sl^ > 0, excitations 
relative to the adiabatic vacuum state can be defined by 
means of correlated quasiparticle annihilation operators, 



.4 



K 



(14) 



which satisfy [j4ki^k , 1 = £kk'- The adiabatic vacuum 
state, |0), is characterized by the absence of quasiparti- 
cles, A K |0) = 0, and is generally a complicated many- 
body state when expressed in terms of the free quanta. 



It is important to recognize that (14) refers to the 
Schrodinger representation. The corresponding Heisen- 
berg operator, A K (t), can be obtained by replacing the 

by the 

(15) 



Schrodinger operators q K and p K entering in 
corresponding Heisenberg operators given in ( 



where A\, is the quasiparticle annihilation operator asso- 
ciated with the initial form of the mass function ^ 2 (x, t\). 
The Bogoliubov coefficients (t) and V£i it) can be ob- 
tained from the basic mode coefficients F R (t) and G A (i) 
by a sequence of elementary transformations. 

Since the operators Ajs.(t) refer to the evolving eigen- 
representation, they are less suitable for the analysis 
which is better performed in the x and k representations, 



A(x, i) = £ [^ L (x, t)i L + V L (x, t)*A\ 

L 

i k (t)=^[z£(x,t)i L + v£(x,t)*4 



(16) 
(17) 



The corresponding Bogoliubov coefficients are obtained 
by the appropriate additional transformations, 



W L (x;t) 



X>|«) S KK (t)Uk(t) 

K.K 

£<k| K ) S KK (t)Uk(t) 

K.K 



(18) 
(19) 



With the above preparations, we are now able to derive 
expressions for any quasiparticle observable. Of particu- 
lar interest is the spatial correlation function, 

p(x,x';i) = (t|it (x )i( x ')| t ) = (t i \A*(x,t)A(x',t)\t i ) 
= ^ u (x,x';i) + p cl (x,x';t) . (20) 

Here the first term arises from the quantum fluctuations 
in the initial state and is given by 



^ u (x,x';t) = ^V L (x,t) V L (x',t)* 



(21) 



The second term is caused by the real agitations in the 
initial state (which are often referred to as the "classical" 
or "thermal" field fluctuations). 

Analogously, we may consider the one-body density 
matrix in the k representation, 



(t\AlA v \t) 



1 k" fl k'l < '/ — V-il-^k 
-qu (±\ | -cl 



Pkk' (*) = 

where the quantal contribution is 



(iil4(*)A'(*)l*i> 



(22) 
(23) 

(24) 



Expressions for observables reflecting two- or many- 
body correlations can also be derived but they become 
progressively more complicated. Moreover, the results of 
course depend on what is being assumed regarding the 
corresponding correlations in the initial state. 
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IV. ILLUSTRATIVE APPLICATION 

In order to demonstrate the practical applicability of 
the described treatment, and at the same time illustrate 
the importance of the quantum fluctuations, we consider 
a one-dimensional system having the most important fea- 
tures expected in actual scenarios. In order to emulate 
the finite transverse extension of systems generated in 
high-energy heavy- ion collisions, we assume that the ini- 
tial mass function exceeds the free value m within only a 
limited region, |x| < 5 fm. We furthermore assume that 
H 2 {x) — to 2 maintains this spatial profile in the course of 
time, while the overall magnitude scales in proportion to 
its bulk value, /Un(t). Fig. [j] shows the time dependence of 
this latter quantity, while Fig. |^ depicts the initial spatial 
profile of the mass function, p(x,0). 



itive throughout (in fact it is never below the free mass 
to), it is possible to analyze the system in terms of the 
adiabatic quasiparticle modes at any time. 

The contribution from the quantum fluctuations to the 
spatial quasiparticle density at various times during the 
evolution, is shown in Fig. |3| and the associated spectral 
distribution is shown in Fig. |], 

P»(x,t) = (P»(x,x;t) = ^,\V L (x,t)\ 2 > ( 25 ) 

L 

Pt(t)=PZ(t) = T,\ V *( t )\ 2 ■ ( 26 ) 

L 

These quantities show what would result if there were no 
excitation at all in the initial state (corresponding to zero 
temperature). Any actual agitations in the initial state 
would then lead to additional contributions. 
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FIG. 1. The time dependence of the bulk value of the 
effective mass function used in the illustration, no(t)- It has 
a rough correspondence with typical mass functions resulting 
from dynamical DCC studies with the linear a model |^J^|. 
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FIG. 3. The spatial quasiparticle density at successive 
times t (as indicated), starting at t — from the correlated 
vacuum associated with the mass function shown in Fig. ^| 
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FIG. 2. The initial mass function used in the illustration, 
(x{x,t = 0). It is obtained by smearing a sharp distribution 
on the interval \x\ < 5 f m with an exponential having a width 
of one fermi. The mass excess p 2 (x,t) — m 2 is assumed to 
retain its spatial profile while the overall magnitude scales as 
the bulk value po (t) given in Fig. |l| 

Given the particular mass function p(x,t), the mode 
functions F x (x,t) and G x (x,t) are then determined nu- 
merically by solving the equation of motion (j^) for each 
canonical mode A. Since the mass function remains pos- 




Momentum ft (MeV/c) 
FIG. 4. The spectral distribution of the quasiparticles at 
successive times (indicated in fm/c) starting at t — from the 
vacuum associated with the mass function shown in Fig. ^| 

At the initial time, t = ti, the quantities p qu (x) and p^ n 
vanish by design. But soon, as the effective mass expe- 
riences its first rapid decrease, energy stored in the high 
initial mass is liberated and a rapidly rising quasiparticle 
density appears with a spatial profile resembling that of 
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the mass excess 5/i(x) = fi(x) — m. Part of these quasi- 
particles are reabsorbed when the mass again increases 
and this cycle is repeated in a less noticable manner, as 
the mass oscillations proceed with an ever decreasing am- 
plitude. The quasiparticles produced by this mechanism 
have a rather structureless spectrum of approximately 
gaussian form. 

Meanwhile, the regularity of the temporal oscillations 
in the mass function begins to manifest itself and, at 
somewhat later times, the spectral profile exhibits a vis- 
ible enhancement of the modes near the resonance fre- 
quency, w res ~ \u a — 300 MeV, corresponding to a mo- 
mentum of k TCS w 266 MeV/c. As a consequence, the 
quasiparticle gas begins to resemble two countermoving 
beams with the corresponding flow velocities centered 
around i>fl ow f=a ±k ros /m. The two corresponding out- 
wards bursts are clearly visible in Fig. |^ at late times 
when the mass function has largely relaxed to its free 
form and the production processes have subsided. 

Traces of these characteristic effects might be observ- 
able in actual collision experiments, such as those under- 
way at RHIC. More extensive calculations would clearly 
be needed for a quantitative assessment of this prospect. 

V. DISCUSSION 



accurate, since the quantum treatment, while having a 
large effect on specific signals, is expected to have a rela- 
tively small overall effect on the effective mass. Neverthe- 
less, an interesting and more formal task, beyond scope of 
present study, is posed by the question of how to best de- 
termine the effective mass function selfconsistently with 
the evolution of the quantum many-body state. 

As noted above, the computational challenge posed by 
the treatment is not prohibitively larger than that met at 
the classical level, since it merely requires (repeated) so- 
lution of the same field equation (|^) . The treatment may 
thus be applied to scenarios similar to those addressed 
at the classical level. This presents an obvious practical 
task which should yield more quantitative information on 
the observable signals of the expected chiral dynamics. 

ACKNOWLEDGEMENTS 

The author is pleased to acknowledge helpful discus- 
sions with D. Boyanovsky, S.J. Jeon, and E. Mottola. 
This work was supported by the Director, Office of 
Energy Research, Office of High Energy and Nuclear 
Physics, Nuclear Physics Division of the U.S. Depart- 
ment of Energy, under Contract DE-AC03-76SF00098. 



In the present paper we have sketched how the non- 
uniform and rapidly evolving scenarios of interest in 
connection with disoriented chiral condensates may be 
addressed with real-time non-equilibrium quantum-field 
theory within the mean-field approximation, in which the 
interactions are encoded into an effective mass function. 
For an arbitrary form of the mass function, the strat- 
egy is to first obtain the time evolution operator, which 
depends only on /x 2 (x, t) , and then evaluate the observ- 
ablcs resulting from any particular initial quantum state 
\t\). This makes it economical to perform averages over 
ensembles of initial states having similar mass functions. 

When the eigenvalues of the mass matrix are posi- 
tive, fl^ > 0, it is particularly instructive to consider 
quasiparticle excitations defined relative to the adiabatic 
eigenbasis. But it should be noted that even if the system 
enters the classically forbidden region where the quasi- 
particles cannot be defined, the dynamical propagation 
of the system works without modification. 

A first application of the developed treatment has been 
made to a system with a non-trivial mass function de- 
pending on both time and space. In addition to demon- 
strating the practicality of the approach, the example 
serves to illustrate the quantitative importance of includ- 
ing the quantum fluctuations in the dynamical treatment. 

The approach assumes that the effective mass function 
is given. For the time being, an approximate form of the 
mass function for a given scenario can be obtained on the 
basis of semiclassical simulations Q. Though not fully 
satisfactory, this approach may actually be reasonably 
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